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1. Introduction 

In recent years, the non-perturbative description of QCD on the lattice has made a significant 
breakthrough in tackling the systematic effects present in the determination of several important 
physical quantities, opening the way for a direct connection to experiments (see e.g. ^ for 
recent reviews). Simulations containing the dynamics of the light-quark flavours in the sea, as 
well as those due to the strange quark and recently also to the charm, using pseudoscalar masses 
below 300 MeV, lattice extents L > 2 fm and lattice spacings smaller than 0.1 fm are presently being 
performed by several lattice groups. Such simulations will eventually allow for an extrapolation of 
the lattice data to the continuum limit and to the physical point while keeping also the finite volume 
effects under control. 

The European Twisted Mass collaboration (ETMC) has carried out large scale simulations with 
Nf = 2 flavours of mass degenerate quarks using Wilson twisted mass fermions at maximal twist. 
Four values of the lattice spacing ranging from 0.1 fm down to 0.051 fm, pseudoscalar masses 
between 280 and 650 MeV as well as several lattice sizes (2.0 — 2.5 fm) are used to address the 
systematic effects. 

The physics of the light pseudoscalar meson is in a suitable sector for investigating the sys- 
tematic effects arising from the continuum, infinite volume and chiral extrapolations of lattice data, 
since the pion mass and decay constant can be measured with high statistical accuracy in lattice 
simulations. Moreover, chiral perturbation theory (#PT) is best understood for those two quantities. 
As a result of this study one can extract several important quantities, such as the u,d quark masses, 
the chiral condensate or the low-energy constants of #PT. First results for the pseudoscalar mass 
raps and decay constant fp$ from these Nf = 2 simulations have been presented in Refs. [§— gl. 

ETMC is currently generating Nf = 2 + 1 + 1 ensembles including in the sea, in addition to the 
mass degenerate light u,d quark flavours, also the heavier strange and charm degrees of freedom. 
First results for the pseudoscalar mass and decay constant from this novel setup were presented 
in [|T§ [n]]. 

In the following we will mainly focus on the results from the Nf = 2 data for raps and fp$. 

2. Lattice Action and Simulation Setup 

In the gauge sector we employ the tree-level Symanzik improved gauge action (tlSym) [Q|. 
The fermionic action for two flavours of maximally twisted, mass degenerate quarks Jl3| , |I4| ] in the 
so-called twisted basis (where the quark field doublets are denoted by % an d %) reads 

S tm = a 4J £{x(x) [Dp] +m + W5^ 3 ] x{x)} , (2.1) 

x 

where rao is the untwisted bare quark mass tuned to its critical value ra cl -i t , \i q is the bare twisted 
quark mass, T 3 is the third Pauli matrix acting in flavour space and D[U] is the Wilson-Dirac 
operator. 

At maximal twist, i.e. rao = ra cr i t , physical observables are automatically O(a) improved with- 
out the need to determine any action or operator specific improvement coefficients [|l4j] (for a 
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Table 1: Ensembles with Nf = 2 dynamical flavours produced by the ETM collaboration. We give the 
ensemble name, the values of the inverse bare coupling j3 = an approximate value of the lattice spacing 
a, the lattice volume V = L 3 • T in lattice units, the approximate value of m^L, the bare quark mass \i q in 
lattice units and an approximate value of the light pseudoscalar mass raps. 



review see Ref. [15]). With this being the main advantage , one drawback of maximally twisted 
mass fermions is that parity and flavour symmetry are broken explicitly at finite values of the lattice 
spacing, which amounts to 0(a 2 ) effects in physical observables. 

For details on the setup, tuning to maximal twist and the analysis methods we refer to Refs. [BL 
0, ^J. Recent results for light quark masses, meson decay constants, the pion form factor, k-% 
scattering, the light baryon spectrum, the v\' meson and the ft) — p mesons mass difference are 
available in Refs. @-||]. 

Flavour breaking effects have been investigated for several quantities ^ ^ ^| 0, |l9| ]. With the 
exception of the splitting between the charged and neutral pion masses, other possible splittings in 
the unitary observables so far investigated are found to be compatible with zero. These results are 
in agreement with a theoretical investigation using the Symanzik effective Lagrangian 24 ] . 

A list of the Nf = 2 ensembles generated by ETMC can be found in table pi. 



3. Scaling to the Continuum Limit 

Here we analyse the scaling to the continuum limit of the pseudoscalar meson decay constant 

'Other properties being, e.g., that the quark mass renormalises only multiplicatively and that the determination of 
the pseudoscalar decay constant does not require a renormalisation factor. 
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Figure 1: (a) ro/ps as a function of (romps) 2 , (b) Continuum limit scaling: ro/ps as a function of (a/ro) 2 at 
two fixed values of romps and at fixed volume. 



fps at fixed reference values of the pseudoscalar meson mass raps and of the lattice size L (we refer 
to [Q, |5|, ||] for details). The purpose of this scaling test is to verify that discretisation effects are 
indeed of 0{a 2 ) as expected for twisted mass fermions at maximal twist. 

In order to compare results coming from different values of the lattice spacing it is convenient 
to measure on the lattice the hadronic scale ro [|25|]. It is defined via the force between static quarks 
at intermediate distance and can be measured to high accuracy in lattice QCD simulations. For 
details on how we measure r^/a and on how its chiral extrapolation is performed, we refer to 
Ref. @i. 

In figure 1(a) we plot the raw data for ro/ps as a function of (roraps) 2 . The vicinity of points 
coming from different lattice spacings along a common curve is an evidence that lattice artifacts 
are small for these quantities. This is indeed confirmed in figure 1(b) where the continuum scaling 
of ro/ps is illustrated: the very mild slope of the lattice data shows that the expected 0(a 2 ) scaling 
violations are small. The result of a linear extrapolation in (a/ro) 2 to the continuum limit is also 
shown. 

Analogous continuum limit scaling studies have been performed for the charged pion mass [Q, 
|9|] and the nucleon mass [19, 26], showing also in these cases signs of only small scaling violations. 



4. ^PT Description of Lattice Data 

The chiral extrapolation of lattice data down to the physical point is currently one of the main 
sources of systematic uncertainties in the lattice results. The possibility to rely on an effective 
theory such as #PT to guide this extrapolation is therefore of great importance to quote accurate 
results from lattice simulations. On the other hand, while smaller quark masses are becoming 
accessible in numerical simulations, the possibility to perform a quantitative test of the effective 
theory as well as to determine the low energy parameters of its Lagrangian becomes more and more 
realistic. 
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We now proceed by presenting the results of a combined continuum, chiral and infinite volume 
extrapolation of mps and /ps for two values of the lattice spacing (corresponding to a 0.063 fm 
and a 0.079 fm). We use the chirally extrapolated values of r^/a to relate data from the two 
lattice spacings and a non-perturbative determination of the renormalisation factor Zp [ |2S| ] in order 
to perform the fit in terms of renormalised quark masses. This analysis closely follows those 
presented in Refs. [Q, ||, to which we refer for more details. 

We perform combined fits to our data for /ps, »ips, ro/a and Zp at the two values of the lattice 
spacing with the formulae: 



roM = rofo 



Kf DH (L) 



;i-2^1og(|) + r / NNLO + D /ps (aAo 
(romps) 2 = [l + % log (jfj + 7r LO + D ms {a/r Q f\ K™(L) 



(4.1) 



v 3 - 

with § = 2fioAi;?/(47r/o) 2 , Xn = 2SoM« , Mr = M<?A>, /o = y/2F Q . T^ denote the continuum 



NNLO terms of the chiral expansion Q29p , which depend on Ai_4 and kjn and kp . The finite size 
corrections factors ^™ H (L) refer to a continuum #PT description using the resummed Liischer 
formula, that we denote as CDH [p0|]. 2 

Based on the form of the Symanzik expansion in the small quark mass region, we parametrise 



in eq. ( |4. 1[ ) the leading cut-off effects by the two coefficients Oy ps mps . Setting Df FSttnps = is equiv- 
alent to perform a constant continuum extrapolation. Similarly, setting T™ LO = corresponds to 
fit to NLO #PT. 

From the fit parameters related to the quark mass dependence predicted by #PT (in particular 
from A3 4, So and /o) the low energy constants £3 4 and the chiral condensate £ can be determined, 
using the following expressions: 



4= log 



A? 



(rnt) 2 ; (42) 

v _ g o/o 2 
L ~ 2 ■ 

By including or excluding data points for the heavier quark masses, it is in principle possible 
to explore the regime of masses in which NLO and/or NNLO SU(2) #PT applies. We have actually 
generalised this procedure in order to estimate all the dominant sources of systematic uncertainties 
that can be addressed from our setup, which include, discretisation effects, the order at which we 
work in #PT or finite size effects. The idea is to use different fit ansatze (see below) on a given 
data-set and to repeat this same procedure over different data-sets: by weighting all these fits by 
their confidence level we construct their distribution and estimate the systematic error from the 
associated 68% confidence interval. 

The fit ansatze we consider are: 

• Fit A: NLO continuum #PT, 7^ LO = 0, D mpsJps = 0, priors for r Ai i2 

• Fit B: NLO continuum #PT, T^ LO = 0, D mpS)/ps fitted, priors for r Ai j2 



For a more detailed description of finite size effects in our data for mps and /ps, we refer to Refs. [ 
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ra M [MeV] 3.54(26) 

h 3.50(31) 

h 4.66(33) 

fo [MeV] 121.5(1.1) 

fx/fo 1.0755(94) 

S [MeV] 2638(200) 

[Zl 1 / 3 [MeV] 270(7) 



Table 2: Summary of fit results. The errors are statistical and systematical, added in quadrature. Bq, £ and 
m u d are (non-perturbatively) renormalised in the MS scheme at the scale ji = 2 GeV. 



• Fit C: NNLO continuum #PT, An PS /ps = 0, priors for r A l 2 and k M F 

• Fit D: NNLO continuum %?T, Om PS ,/ps fitted, priors for rrjAi^ and k.M,F 

The choice of the different data-sets (each of them including data for different lattice spacings, 
quark masses and physical volumes) is made in order to quantify how the quality of the fit is 
modified when including/excluding data from e.g., a given mass region or with a given volume. 
The data-sets considered in the fits are listed in Ref. [g]. 

The parameters &m,f and roAi^ appearing at NNLO (and in the latter case, also at higher 
order in the CDH expressions) need to be fitted with some additional input information (priors). 
Additional lattice data would be needed to allow these parameters to remain free in the fit. The 



parameters rnAi^ use priors from the estimates for £12 in Ref. []30|], while very mild priors are 
used for Um^f = ± 10. For a detailed description of the use of priors and, in general, of the 
statistical analysis, we refer to Ref. 

The results of the combined fits to data from the two lattice spacings a k, 0.063 fm and 
a m 0.079 fm are given in table [| The errors are statistical and systematical, added in quadrature. 
For a review on the determination of the low-energy constants and of the light-quark mass we refer 
to |32[ , ||]. Recent determinations from ETMC of the pseudoscalar decay constant and chiral 



condensate in the £ -regime were presented in [|33fl. 



As a further check, we have performed additional fits including either the finer lattice spacing 
(a ps 0.051 fm) or the coarser one (a O.lOOfm) finding total compatibility with the results of 
table g 



5. Discussion and Conclusion 



Here we collect a short list of observations coming from a set of #PT fits. A complete descrip- 
tion of these fits was presented in Ref. 

We observe that including in the fits pseudoscalar masses nips > 520 MeV decreases signifi- 
cantly the quality of the NLO fits (^ 2 /dof > 1). This indicates that the applicability of NLO #PT 
in that regime of masses is disfavoured. 

On the contrary, extending the fit-range to a value of raps ~ 280 MeV preserves the good 
quality of the fit and gives compatible values for the fit parameters. This result makes us confident 
that the extrapolation to the physical point is trustworthy. 



6 



Scaling and %PT Description of Pions from Nf = 2 tmQCD 



Gregorio Herdoiza 



Including lattice artifacts in the fits gives results which are compatible to those where D m?s j^ s 
is set to zero but with a somehow better # 2 /dof. We observe that the values of the fitting parameters 
D mps j ps are compatible with zero within two standard deviations. This is in line with the small 
discretisation effects observed in the scaling test. 

The inclusion of NNLO terms produces similar results to the NLO fits in the quark mass 
region corresponding to wips G [280,520] MeV. When fitting data only in this mass region (i.e. 
when excluding from the fit the heavier masses at nips ~ 650 MeV), we observe that the fit curve 
at NLO lies closer to those data points (heavier masses) than the NNLO one. On the other hand, 
when including the heavier masses, the NNLO fit is able describe these data points but the quality 
of the fit is somehow reduced with respect to the NLO fit. To improve the sensitivity of our lattice 
data to #PT at NNLO, additional data points would be needed. 

We have presented determinations of fps and mps and of their continuum, infinite volume 
and chiral extrapolations. As a result, we obtain accurate determinations of the u,d quark mass, the 
chiral condensate as well as of low-energy constants of the effective theory, including an exhaustive 
estimate for the systematic uncertainties. 

The only systematic effect which cannot be addressed by this study corresponds to effect of the 
strange and the charm quarks in the sea. We are currently in the process of extending this analysis 



to the Nf = 2 + 1 + 1 ensembles which are being generated by ETMC [10, |1 lfl. 



We thank the members of ETMC for the most enjoyable collaboration. 
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